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ABSTRACT 

We report on deep imaging of a remote M31 globular cluster, MGCl, obtained with 
Gemini/GMOS. Our colour-magnitude diagram for this object extends 5 magni- 
tudes below the tip of the red giant branch and exhibits features consistent with an an- 
cient metal-poor stellar population, including a long, well-populated horizontal branch. 
The red giant branch locus suggests MGCl has a metal abundance [M/H] « —2.3. We 
measure the distance to MGCl and find that it lies ~ 160 kpc in front of M31 with 
a distance modulus fi = 23.95 ± 0.06. Combined with its large projected separation 
of Rp = 117 kpc from M31 this implies a deprojected radius of i?gc = 200 ± 20 kpc, 
rendering it the most isolated known globular cluster in the Local Group by some 
considerable margin. We construct a radial brightness profile for MGCl and show 
that it is both centrally compact and rather luminous, with Afy = —9.2. Remarkably, 
the cluster profile shows no evidence for a tidal limit and we are able to trace it to 
a radius of at least 450 pc, and possibly as far as ^ 900 pc. The profile exhibits a 
power-law fall-off with exponent 7 — —2.5, breaking to 7 = —3.5 in its outermost 
parts. This core-halo structure is broadly consistent with expectations derived from 
numerical models, and suggests that MGCl has spent many gigayears in isolation. 

Key words: globular clusters: general, galaxies: individual: M31 



1 INTRODUCTION 

Cosmological models of structure formation predict that 
galaxies are built up via the hierarchical accretion and 
merger of many smaller subsystems over time. It is expected 
that the signatures of these galaxy assembly processes may 
be seen in the outskirts of present-day large galaxies where 
dynamical timescales are very long. Globular clusters, as lu- 
minous compact objects that are found out to distant radii 
in the halos of massive galaxies, are thus potentially ex- 
tremely useful probes of galaxy formation. This is especially 
the case for systems beyond a few Mpc where globular clus- 
ters are often the only accessible tracers of ancient stellar 
populations. 

In terms of testing and refining these ideas, galaxies 
in the Local Group play a central role. They are close 
enough that their various components, including individ- 



ual globular clusters, may be studied in detail as resolved 
stellar populations. M31 is of particular significance as it 
is the nearest large galaxy to our own, and possesses the 
most extensive globular cluster system in the Local Group. 
This galaxy has been targeted by a number of wide-field 
imaging surveys in recent years and these have revealed a 
wealt h of low-brightnes s substructure in its outer regions 



(e.g.. Ilbata et all I2OOII: Ferguson et"all I2OO2I : llbata et all 
l2007l : iMcConnachie et al.ll2009[ ) as well as many new mem- 
bers of its halo globular cluster system, extending to very 
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Although the M31 globular cluster system shares many 
similarities with that of the Milky Way, as the sample of 
outer halo M31 globular clusters has grown and been stud- 
ied in more detail over the past few years several clear dif- 
ferences have emerged. Among these are the existence in 
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the M31 system of very extended yet still rather luminous 
clust ers l|Huxor et al ] |2005l : lMackev et al ] |2006l : lHuxo~ et al.l 
that do not seem to have matching counterparts in the 
Milky Way, and indications that M31 may possess many 
more very luminous compact globular clusters in its outer 
regions than does the Milky Way (e.g.. iMackev et "al ] |2007l : 
iHuxor et al]|2009l ). The origin of these distinctions is not 
yet clear but it has been speculated that they may reflect 
differences in the formation and/or accretion histories of the 
Milky Way and M31. It is therefore of considerable interest 
to obtain as much additional information as possible about 
the remote globular cluster system in M31. 

One of the most intriguing newly- found objects in 
the M31 halo is a b right globular cluster discovered 
by iMartin et al.l l|2006l ) in survey observations from the 
MegaCam imager on the Canada-France-Hawaii Telescope 
(CFHT). The cluster lies at a = 00:50:42.45, 5 = 
-1-32:54:58.70 (J2000.0), or roughly ~ 8.5° from the centre 
of M sfl Adopting the standa rd distance to M31 of « 780 
kpc (|McConnachie et al.ll2005h . this angular separation cor- 
responds to a projected radius of Rp = 117 kpc, placing the 
new globular cluster among most remote such objects yet 
known in the M31 system and, indeed, in the Local Group 
as a whole. Following the nomenclature adopted in the lat- 
est version (V3.5) of the Revised Bologna Ca talogue (RBC) 
of M31 globular clusters (jCalleti et al.l[2007[ ). we hereafter 
refe r to this objec t asMGCl. 

iMartin et all (j2006h used their CFHT/MegaCam dis- 
covery imagery to measure a few characteristics of the 
cluster. They found it to be a compact, luminous object 
with half-light radius ~ 2.3 ± 0.2 pc and absolute inte- 
grated magnitude My = -8.5 ± 0.3. Their CMD showed a 
well-populated red giant branch, and a probable horizontal 
branch near their faint detection limit. Isochrone fltting im- 
plied a metal abundance [Fe/H] ~ —1.3. Most intriguingly, 
the observed brightness of the tip of the red giant branch, 
along with that of the putative horizontal branch, suggested 
that MGCl may lie considerably closer than M31, meaning 
that the cluster may be extremely isolated indeed with a 
deprojected galactocentric radius of 175 ± 55 kpc. 

Two gro u ps ha ve published radial velocities for MGCl: 
iGalleti et al] |20o3) obtained Vr = -312 ± 17 kms~^ frorn 
low-resolution spectroscopy, while lAlves-Brito et all (|2009l ') 
measured Vr = —354.7 ± 2.2 kms~^ from a high resolution 
Keck spectrum. These measurements are consistent with the 
systemic velocity of M31 {vr — —300 kms^^), suggesting 
that MGCl is a bona fide member of M31's outer halo glob- 
ular cluster system. In this regard it is of considerable sig- 
nificance as very few such objects have yet been discovered 
- there are only thr ee other known globular clusters in M3 1 
with Rp > 70 kpc (|Mackev et al.ll2007l : iHuxor et al]|2008h . 
Indeed there are also very few individual halo stars that 
have been iden tified at projected radii comparable to t hat 
of MGCl (e.g., lKalirai et al.ll2006l : IChapman et al.l l2008l. 

In this paper we present results from deep Gemini 
follow-up imaging of this intriguing cluster, with the partic- 
ular goals of (i) obtaining an accurate distance measurement 
in order to understand properly its position in the M31 sys- 



1 Assumed to be at a = 00:42:44.31, <5 = -1-41:16:09.4 (J2000.0). 



tem; (ii) obtaining a precise estimate of its metal abundance; 
and (iii) making a thorough analysis of its structure. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Imaging and photometry 

Our deep imaging of MGCl was carried out in queue mode 
on the night of 2007 October 13 using the Gemini Multi- 
Object Spectrograph (CMOS) at the 8.1m Gemini North 
telescope on Mauna Kea, Hawaii. The data were obtained 
in clear, photometric conditions and under excellen t seeing 
(0.4 - 0.5")- The CMOS imager l|Hook et al.ll2004h is com- 
prised of three adjacent 2048 x 4096 pixel CCDs separated 
by gaps of ~ 2.8", and has a field of view (which does not 
cover the full CCD package) of 5.5' x 5.5'. To take advan- 
tage of the high quality conditions we employed unbinned 
imaging, resulting in a plate scale of 0.0727 arcsec/pbcel. We 
obtained our observations in the CMOS g' , r' , and i' filters, 
which are very similar to the g, r, and i filters used by the 
Sloan Digital Sky Survey (SDSS ). Details of the SDSS pho- 
tometric system may be found in lFukugita et akl (|l996l ). We 
took six images per filter, arranged in a 3 x 2 dither pat- 
tern with a step size of 5" designed to eliminate the gaps 
between the CCDs and provide continuous coverage of our 
field. Exposure durations were 400s per image for the r' and 
i' filters and 600s per image for the g' filter. 

We reduced our data using the CMOS reduction pack- 
age in IRAF. Appropriate bias and flat-field images obtained 
as part of the standard CMOS baseline calibrations were 
downloaded from the Gemini science archive and then ap- 
plied to each exposure using the GIREDUCE task. The three 
CCD frames in a given exposure were next mosaicked into a 
single frame using the GMOSAIC task, and the six mosaicked 
frames for a given fllter were then combined into a single 
image using the imcoadd task. 

Fig. [T] shows the central 70" x 70" region of our com- 
bined 2400s i'-band image. Stellar objects in this image have 
a FWHM ^ 6.2 pixels, or 0.45", and the outer region of 
MGCl is clearly resolved into stars. The inner part of the 
cluster is unresolved, but, as described below, the images 
nonetheless contain a sufficient number of individual clus- 
ter members to produce a well-populated colour-magnitude 
diagram (CMD). 

We performed photometric measurements on our com- 
bined images using the stand-alone versions of daophot 
11 and ALLSTAR II (Stetson 1987). Because of the crowded 
nature of the region of interest, full point-spread function 
(PSF) fitting photometry was necessary. For each combined 
image we used daophot ii to conduct a first pass of object 
detection, and selected ~ 75 relatively bright, isolated stars 
to construct an initial PSF. Some experimentation showed 
that the best fit was obtained by selecting a Gaussian PSF 
model and allowing quadratic variation across the field of 
view. Stars with an error in their fit more than three times 
the average were removed from the list, and the PSF re- 
defined. After iterating until convergence, at which point 
there were typically still ~ 50 — 60 stars defining the PSF 
model, we used ALLSTAR to subtract from the image all stars 
in close proximity to those used to construct the PSF. The 
now-isolated PSF stars on this subtracted image were used 
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to recalculate and further refine the PSF model. Next, we 
used ALLSTAR to apply this model to the original image and 
subtract all known stars. This subtracted image was then 
run through daophot in order to find faint objects missed 
in the first detection pass. We then took the original im- 
age and our final PSF model and used ALLSTAR to perform 
photometric measurements on the complete list of detected 
objects. 

In order to eliminate non-stellar objects or objects with 
poor photometry, we passed the resulting list of measure- 
ments through several quality filters based on information 
calculated during the PSF fitting. Specifically, we filtered 
objects by the their fit, their estimated photomet- 

ric error, and their measured sharpness relative to the PSF 
model. To define suitable ranges in these parameters for well- 
measured point sources we selected all objects in an annulus 
of radius 12" — 40" centered on the cluster, corresponding to 
the region populated almost exclusively by cleanly resolved 
member stars. These objects defined clear, narrow loci in 
terms of x^, photometric error and sharpness as functions 
of magnitude in each of the three filters, which we used to 
filter the final list of detections across the full field of view. 

2.2 Photometric calibration 

In addition to observing our target field we also obtained 
imaging of photometric stan dard stars in one equatorial 
field, around star SA 92-249 (iLandoltHigg^ l. This imaging 
was conducted through all three filters immediately after the 
observations of our primary field were complete. We took 
two exposures per filter, with durations of 4s, 3s, and 2s 
each in g' , r', and i' , respectively. An identical reduction 
process to that followed for the primary imaging was applied 
to the observations of our photometric standards to produce 



one mosaicked, combined image per filter. We then used the 
stand-alone version of daophot ii to perform aperture pho- 
tometry on the standard stars in these combined images, 
using an aperture equal in radius to that for the PSF mod- 
els we used to obtain photometric measurements from the 
main cluster images. 

In all, six photometric standard stars were covered by 
the GMOS field of view. Unfortunately, however, we could 
only obtain accurate photometric measurements for four of 
these - one star was completely vignetted by the GMOS On- 
Instrument Wavefront Sensor (OIWFS), which facilitates 
precise guiding and use of the telescope's active optics sys- 
tem, while a second star fell on the gap between the CCDs 
in half the images. 

Our standard field lies within the SDSS footprint; how- 
ever the stars are sufficiently bright that they are fiagged 
as saturated in the SDSS catalogue. We therefore employed 
an alternative method of obtaining their S DSS r-, and 
i-band AB-magnitudes, by transforming the lLandoltl (|l992l ) 
Johnson-Kron-Cousins BVRI measurements using the re- 
lations listed on the SDSS Data Release 7 (DR7) websitfl 
Specifically, we used inverted versions of Lupton's ugriz to 
BVRI transformations: 

g = 0.6489 B-f 0.3511/ - 0.1460 

r = 1.3088 ii- 0.3088/ 4- 0.0704 (1) 
i = 0.2570 /? + 0.7430 / -I- 0.3208 

The other alternative transformations listed 
on the DR7 w ebpage, in particular those of 
IJester et all (120051 ): iBihr. Karaali fc Tuncell (120051 ): and 



\http://www.sdss. org / dr7 / algorithms /sdss UB VRI Trans form.html\ 
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Figure 2. Transformations from GMOS instrumental magni- 
tudes to SDSS AB magnitudes derived from the observed pho- 
tometric standards. The four stars with complete measurements 
are marked as solid points, while the star which fell on the CCD 
gap in half the images is an open point. 

Ijordi. Grebel fc Ammoril (|2006h give matching gri mag- 
nitudes for our standards to within a few x 0.01 mag. 
Therefore, the selection of any of these sets of relations 
would not significantly alter our derived photometric 
cafibration. 

Fig. [2] shows the relationship between our measured 
GMOS instrumental g' , r' and i' magnitudes for the stan- 
dard stars, and their derived SDSS gri AB-magnitudes. The 
instrumental magnitudes as plotted have been corrected to 
the exposure times of the three respective combined cluster 
images. In each panel the star marked with an open circle 
was that which fell on the gap between the CCDs in half 
the images. We did not use our measurements for this star 
when calculating the best-fit transformations. From Fig. [2] 
the relationships between g'r'i' and gri are well described 
by a simple zero-point shift. We determined the best-fit cal- 
ibration via linear regression: 

g = 3' + 9.84 

r = r' -I- 9.75 (2) 
i = i' + 9.50 

Formal (random) uncertainties on these zero-point shifts 
from the linear regression are typically ±0.02 mag. Even 
so, the above transformations were derived using only four 
stars which do not span a particularly large range in colour 
(0.45 ^ (5 — j) ^ 1.35), meaning that it was imposs i ble for 
us to determine colour terms for Eq. [S] Ijorgenseni (|2009l ) 
has demonstrated that small colour terms are expected in 
transformations to the SDSS photometric systems for all 
GMOS-N filters, and has used a large number of standard 
star observations to derive these corrections. Unfortunately 
her photometric calibration and colour terms are valid only 
for data obtained prior to November 2004 when the main 



mirror of Gemini North was re-coated with silver in place 
of aluminium, and so we cannot use them directly for our 
MGCl observations. None the less, we expect them to give 
some indication of the size of the systematic errors intro- 
duced into our photometric calibration by not including any 
colour terms in Eq. [5] 

Our colour-magnitude diagrams show that stars in 
MGCl span -0.2 < {g - r) < 1.4 and -0.4 < {g - i) < 1.7 
(see Section I3.1|l . The appropriate colour-term corrections 
to gri magnitudes from I JorgensenI (j2009l ) are then: 
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where the two corrections Ami and Am2 for a given filter 
may be calculated independently. Over the range of colours 
spanned by cluster members, these would correspond to sys- 
tematic offsets in our SDSS magnitudes of —0.04 < A5 < 
0.04, -0.02 < Ar < 0.02, and -0.03 < Ai < 0.03, and 
in our SDSS colours of -0.02 < A{g - r) < 0.02 and 
—0.01 < A(5 — i) < 0.01. Such uncertainties are comparable 
in size to the formal random uncertainties on our zero-points 
(Eq. [2]). As such, we are confident that in what follows our 
analysis and conclusions are largely unaffected by the fact 
that we did not include colour terms in our photometric cal- 
ibration. For completeness, where appropriate we explicitly 
make note of the effects these small systematic offsets may 
have on our measurements of MGCl. 



3 ANALYSIS OF THE CMD 
3.1 General properties 

Our final (5, r) and (5, i) colour-magnitude diagrams are pre- 
sented in Fig. [3] The upper panels show all detected objects 
which passed successfully through the quality filter we de- 
scribed in Section 12.11 while the lower panels show the sub- 
set of those objects which lie in the projected radial range 
12" - 40" from the centre of MGCl. We determined these 
radial limits empirically, to exhibit the best-defined clus- 
ter sequences with minimal degradation due to crowding or 
background contamination. The lower panels also show the 
typical photometric uncertainties as functions of magnitude. 

The cluster sequences are clearly visible in the upper 
panels in Fig. [3] as are several contaminating populations. 
To the red side of the CMDs, at (5 - r) > 1.2 and (5 - i) > 
1.6, lies the sequence due to Galactic disk dwarfs; this is 
most prominent in the (g-,r) CMD. Also evident are a few 
members of the Galactic halo (e.g., at {g — i) 1.0 and i < 
22.5) and a significant number of faint unresolved galaxies 
(which have 0.0 < (5 - «) < 1.4 and i > 25.0). In the lower 
panels, which cover only a limited area on the sky, objects 
from these populations are almost completely absent. 

The cluster itself possesses a well-populated steep red- 
giant branch (RGB) and a horizontal branch (HB) that ex- 
tends far to the blue. Our photometry reaches ~ 1.5 mag- 
nitudes below the HB, meaning that this feature is clearly 
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Figure 3. Final {g, r) and {g, i) colour-magnitude diagrams for the full GMOS frame (upper panels) and for an annulus spanning the 
radial range 12" — 40" centred on MGCl (lower panels). Typical photometric uncertainties are marked in the lower panels. 



isolated on the CMD. There are no large gaps along the mea- 
sured length of the HB, implying that the cluster is likely to 
possess a population of RR Lyrae stars. Our observational 
baseline, however, does not span a sufficiently long interval 
to allow a useful variability search to be performed. Overall, 
the morphology of the HB and the implied presence of RR 
Lyrae stars, together with the steep RGB, is indicative of a 
very old (> 10 Gyr) metal-poor stellar population. 

For ancient systems such as this, the locus of the RGB 
above the level of the HB is mostly dependent on the metal 
abundance and only weakly sensitive to age. The observed 
shape of this feature therefore provides a usefu l and well- 
estab lished probe of the cluster metallicity (e.g.. [Saraiedinil 
1 19941 ). Meanwhile, the observed level of the HB and the ap- 
parent colour of the RGB at the HB level allow estimates 
to be made of the distance to the cluster and the amount of 
foreground reddening. In order to obtain such measurements 
for MGCl, we compared our CMDs with fiducial sequences 
observed for several Galactic globular clusters in the SDSS 
photometric system. We subsequently checked the consis- 
tency of our results by fitting theoretical isochrones from 
stellar models calculated by two different groups. 



3.2 Comparison with Galactic globular clusters 

lAn et all (|2008l ') provide CMDs in the SDSS photometric 
system for 17 Galactic globular clusters spanning a wide 
range of metal abundances, —2.3 < [Fe/H] < —0.7. Unfortu- 
nately for the purposes of estimating a photometric metallic- 



ity for MGCl, many of the CMDs they provide (and in par- 
ticular those for all but one of the clusters with [Fe/H] < —2) 
are truncated only a magnitude or two above the HB level 
due to the saturation limits of the SDSS imagery (see their 
Fig. 12). None the less, there are two well-populated clus- 
ters with abundances in the range of interest that are suf- 
ficiently distant for the entire extent of their RGBs to be 
sampled. These are NGC 2419 and NGC 7006, whic h have 
[Fe/H] = —2.12 and —1.63, respectively, according to iHarrij 
(jl999 ). 

To correctly align the fiducial sequences for these two 
clusters, we shifted them vertically on the CMD so that 
their HB levels matched that observed for MGCl, and hor- 
izontally so that their RGB colours at the level of the HB 
matched the observed value for MGCl. We determined the 
appropriate HB level for MGCl in each of the three pass- 
bands by plotting luminosity functions. These are presented 
in Fig. 13] - the HB level is clearly visible in each of the 
histograms, at g = 25.1, r = 24.9, and i = 24.7. The un- 
certai nty in these values is ±0.05 mag. Note that I An et all 
(|2008t ) do not provide fiducial HB sequences for their clus- 
ters; however their photometry is publicly available and 
these ridgelines are straightforward to determine. 

Because lAn et ahl (|200^) obtained their fiducial se- 
quences from a homogeneous data set via a uniform reduc- 
tion and photometry pipeline, we were able to impose an ad- 
ditional constraint in our alignment procedure - that the im- 
plied offsets in distance modulus (jj,) and foreground redden- 
ing between MGCl and a given reference cluster should be 
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Figure 4. Cluster luminosity functions in each passband, con- 
structed using only stars in the radial range 12" — 40". The HB 
level is clearly visible and is marked with an arrow in each panel. 



equal on both the [g, i) and {g, r) CMDs. That is, for a given 
reference cluster we tested all likely combinations of A/x and 
AE{B — V), converted these to vertical and horizontal off- 
sets in the observational plane, and determined which com- 
bination gave the closest alignment between the observed 
globular cluster fiducial and MGCf on both the {g,i) and 
{g, r) CMDs simultaneously. To convert A/i and AE{B — V) 
to vertical and h o rizont al shifts we used the coefficients from 
IStoughton et all (|2002l ). which relate the foreground extinc- 
tion to E{B - V) in SDSS passbands: Ag = 3.793E{B - V), 
Ar = 2.7^1E{B-V) and Ai = 2m&E{B-V). Given these, 
it is trivial to show that E{g - «) = l.7mE{B - V) and 
E{g -r) = f.042£(B - V). 

The results for NGC 24f 9 and NGC 7006 are presented 
in Fig. [5] The RGB fiducial for NGC 7006 is clearly much 
too red, implying that MGCl is nowhere near as metal rich 
as this cluster. In contrast the fiducial for NGC 2419 pro- 
vides a much better fit along the majority of the RGB, only 
diverging marginally to the red towards the top. This sug- 
gests that MGCl has a metal abundance close to, but still 
below, that of NGC 2419 at [Fe/H] = -2.12. 

To explore t his further, we used the observations of 
IClem et al.l (I2008D who provide a fiducial RGB sequence for 
the very metal-poor Galactic globular cluster M92 ([Fe/H] — 
—2.28) in the u'g'r'i'z' photometric system, which is very 
similar (but not identical) to the ins t rumen tal SDSS ugriz 

I 2008h . The following 



system used here and bv lAn et al 



relations, pro vided by iTucker et al 



( 20061 ). allowed us to 



transform the lClem et~al] (120081 ) M92 fiducial into the SDSS 
ugriz system: 



g = g' + 0.060 [{g' - r') - 0.53] 
r = r' + 0.035 [(r' - i') - 0.21] 
i = i' 0.041 [(r' - i) - 0.21] . 



(4) 



lAn et ai] (|2008l ) demonstrate that these transformations are 



very accurate (better than 2 per cent) for M92 and other 
Galactic globular clusters. Unfortunatelv. lciem et al] (|2008l ) 
do not provide HB fiducials for their clusters, or a public 
archive of their photometry. Therefore, we had to use t he HB 
fiducial for M92 that we derived from the lAn et al.1 (|2008l ) 
photometry - this is evidently not ideal, but the demon- 
strated accuracy of the transformations listed in Eq.|4]means 
that any systematic uncertainties introduced should be, at 
most, a few hundredths of a magnitude. 

The results we obtained by aligning these M92 fidu- 
cials are also presented in Fig. [S] On both the {g, i) and 
((?, r) CMDs the agreement is very close along the full length 
of the RGB for identical distance modulus and foreground 
reddening. It thus seems likely that MGCl has a very low 
metal abundance, with [Fe/H] ~ —2.3. It is worth noting 
th at this value is more metal poor than the initial estimate 
bv iMartin et all (HoOg), who obtained [Fe/H] ~ -1.3, and 
also the value of [Fe/H] ~ —1.37 ±0.15 measured fro m high- 
resolution spectroscopy bv lAlves-Brito et al.1 l|2009l ). 

In Section [2.21 we noted the lack of colour terms in our 
photometric calibration. It is important to assess what im- 
pact this may have on the shape of the RGB and hence on 
our derived metal abundance. On the {g,i) CMD the part 
of the cluster RGB above the HB spans the colour range 
1-0 (p — *) I-'''. According to Eq. [3]the expected sys- 
tematic offsets due to the omission of colour terms span the 
range 0.0 < A{g — i) < 0.01 over this colour interval. Simi- 
larly, the expected systematic offsets on the {g, r) CMD span 
the range 0.01 < A{g — r) < 0.02. In both cases the result of 
adding colour terms is to stretch the RGB to the red, mak- 
ing it appear more metal-rich than in the present CMDs. 
However, the expected magnitude of this effect is too small 
to significantly alter our derived metal abundance. 

The horizontal and vertical shifts required to align 
the M92 fiducials with our CMDs allow us to estimate 
the fo reground redd ening and distance modulu s for MGCl. 
Both iHarrid l|l996t ) and ISchlegel et all l| 19981 ) agree that 
E{B - V) = 0.02 for M92, resulting in an extinction- 
corrected distance modulus of /i = 14.58 for this cluster 
(HaiTis 199^ ). The horizontal and vertical offsets we derived, 
as noted in Fig. [S] then imply that the foreground redden- 
ing towards MGCl is E{B - V) = 0.17, and the distance 
modulus for this cluster is p = 23.93. 

There ar e three additiona l Galactic globular clusters in 
the sample of I An et al.l (|2008l ) that have metal abundances 
comparable to that of M92: M15, NGC 5053, and NGC 5466. 
These all have truncated RGB sequences and are hence not 
useful for obtaining a precise metallicity estimate for MGCl; 
however their observed fiducials are none the less suitable for 
determining additional estimates of the foreground redden- 
ing and distance modulus of MGCl, on the basis that they 
are all sufficiently similar in metal abundance for our align- 
ment procedure to provide physically meaningful results. 

Table [T] summarizes the values of distance modulus and 
foreground extinction derived from fitting these cluster fidu- 
cials to our CMDs; we have also included NGC 2419 for com- 
parison. Random uncertainties in individual values of /j, and 
E{B — V) are related to the errors associated with aligning 
the fiducials correctly. Experimentation shows that fiducials 
shifted by ±0.02 mag in colour and in the range ±0.05 — 0.1 
mag in brightness result in noticeably poorer fits. The rela- 
tionships between E{B — V), and E{g~i) and E{g — r) mean 
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Figure 5. Fiducial sequences (solid lines) for the Galactic globular clusters NGC 7006 (left), NGC 2419 (centre), and M92 (right) 
overlaid on the (g, i) and (g, r) CMDs for MGCl. The RGB sequence for NGC 7006 ([Fe/H] = -1.63) is much too red to provide a good 
fit to that of MGCl. In contrast, the sequence for NGC 2419 ([Fe/H] = —2.12) provides a much better fit but diverges slightly to the 
red on the uppermost RGB, while the sequence for M92 ([Fe/H] = —2.28) provides an excellent fit along the full extent of the RGB. 



Table 1. Results from fitting metal-poor Galactic globular cluster fiducials to the CMD of MGCl. 



Cluster 


[Fe/H]" 


E(B - Vf 




/\E{B - V) 




E{B - 1^)mgci 


MMGCl 


NGC 2419 


-2.12 


0.09 


19.69 


0.09 


4.28 


0.18 


23.97 


NGC 5466 


-2.22 


0.01 


15.97 


0.14 


7.97 


0.15 


23.94 


M15 


-2.26 


0.10 


15.06 


0.09 


8.80 


0.19 


23.86 


M92 


-2.28 


0.02 


14.58 


0.15 


9.35 


0.17 


23.93 


NGC 5053 


-2.29 


0.03 


16.10 


0.16 


7.75 


0.19 


23.85 



° Metal abundance as listed bv lHarrid l|l996l 'l in the 2003 online update^ 

^ Foreground reddening is calculated as the average of the values from iHarrisI lll996l ) and ISchlegel et al.l l|l998 h ; extinction-corrected 
distance modulus is calculated using this value and the apparent visual distance modulus listed bv lHarrisI Jl996V 



that the random uncertainties in E{B — V) are then better 
than 0.02 mag. The contribution to the uncertainty in /i 
from the vertical ahgnment precision is therefore dominant 
over the contribution from the uncertainty in the calculated 
extinction, and we hence assume the random uncertainty in 
distance modulus to be ±0.1 mag. We note that there is 
also a systematic contribution due to the uncertainty in our 
estimation of the correct HB levels - Fig. |4] suggests this is 
of order ±0.05 mag at most. 

Our derived colour excesses and distance moduli for 
MGCl both exhibit small dispersions across the sample 
of comparison clusters. The rms values are E{B ~ V) = 



0.18 ± 0.02 and ^l = 23.91 ± 0.06. We discuss these results 
and their implications in more detail in Section [3.41 below. 

3.3 Comparison with theoretical isochrones 

To check the consistency of our results, including their 
robustness to variations in both the assumed age and 
Q-element enhancement of MGCl (neither of which are 
tightly constrained by observation), we employed theoreti- 
cal isochrones calculated by two different groups. We down- 
loaded models of both scaled-solar and a-enhanced compo- 
sition - Galactic globular clusters typically exhibit enhance- 
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Figure 6. Red-giant branches from scaled-solar isochrones calculated by tlie BaSTI and Dartmouth groups overlaid on our cluster (g, i) 
and {(/, r) CMDs. In each panel the isochrone with the best-fitting metallicity is plotted with a thick red line. The black dotted lines 
mark the adopted cluster HB level, from Fig. |4] 



ment in a-elements compared to the s olar abundance, of 
roue hly [a/Fe] ~ +0.3 (see e.g.,|c arne 

3[l996i)- We initially 
selected models as close as possible to ~ 12.5 Gyr old (taken 
to be representative of the age of the oldest metal-poor 
Galactic globular clusters) , but subsequently tested how our 
results varied with changes in age of ±2 Gyr either side of 
this value. In more detail, we used the following: 

• BaSTI models (|Pietrinferni et al.|[20oi . l2006l ). for which 
we downloaded isochrones of age 12.5 Gyr, scaled-solar com- 
position, and metal abundances of [Fe/H] = —1.79 and 
—2.27; along with isochrones again of age 12.5 Gyr, a- 
enhanced composition with [a/Fe] = -1-0.4, and iron abun- 
dances of [Fe/H] = -2.14 and -2.62. 

• Dartmouth models (jPotter et al.ll2007l . [20o3 '). for which 
we downloaded isochrones of age 12.7 Gyr, metal abun- 
dances sampling the range —2.5 ^ [F^/H] ^ —1.6 at 0.3 
dex intervals, canonical Helium abundance {Y « 0.245), and 
both scaled-solar and a-enhanced composition - the latter 
with [a/Fe] = -1-0.4. Unlike the BaSTI isochrones which in- 
clude post-RGB evolution, the Dartmouth group provides 
separate synthetic HB models. 

Although models by the Padova group l|Girardi et al.ll200d . 
l2004l : llVlarigo et ahlBoO^ ) are also commonly employed, we 
elected not to use these in the present work for two rea- 
sons. First, Padova models with a-enhanced composition are 



not presently available for the suspected age and metallic- 
ity of MGCl; and second, it is known that the scaled-solar 
Padova models in the SDSS ugriz system for ancient metal- 
poor systems generally possess redder RGB loci than should 
be expected. A par ticular example of this may be seen in 
iGirardi et al.l ([2004) where Padova isochrones in the SDSS 
photometric system with [Fe/H] = —1.7 provide a good fit to 
the RGB of the Galactic globular cluster NGC 2419, which 
is ~ 0.4 dex more metal-poor than this. 

We overlaid the isochrones using the same methodol- 
ogy as that described in Section 13.21 - shifting them ver- 
tically on the CMD so that their HB levels matched that 
observed for MGCl, and horizontally so that their RGB 
colours at the level of the HB matched the observed value. 
In this case however, since we were dealing with theoretical 
isochrones (which represent the output of stellar evolution 
models transformed onto the observational plane), we were 
forced to relax the additional constraint that the implied 
distance modulus and foreground extinction be identical for 
both the (gr,i) and (g,r) CMDs. 

Fig.[6]shows aligned scaled-solar isochrones with a vari- 
ety of metallicities on our {g,i) and [g,r) CMDs. The results 
are in excellent agreement with the metallicity estimate of 
[Fe/H] « —2.3 that we obtained by aligning Galactic globu- 
lar cluster fiducials on the CMD. Here, the BaSTI isochrone 
with [Fe/H] = -2.27 fits the cluster RGB most closely, 
whereas the next available isochrone with [Fe/H] = —1.79 is 




clearly too metal rich. Similarly, the Dartmouth isochrones 
with [Fe/H] = -2.5 and [Fe/H] = -2.2 both provide reason- 
ably good matches to the observed RGB, while the isochrone 
with [Fe/H] = —1.9 is clearly too metal rich. From these 
models we therefore infer [Fe/H] ~ —2.35. To confirm that 
such measurements are not sensitive to the selected age we 
repeated the alignment procedure using sets of 10 and 15 
Gyr old isochrones, obtaining consistent results to within 
« 0.15 dex. 

Fig. [7] shows aligned a-enhanced isochrones on our 
((?, i) and (g, r) CMDs. The degree of a-enhancement is 
marginally higher than the [a/Fe] ~ +0.3 typically observed 
for Galactic globular clusters and thus likely represents the 
maximum effect Q-enhancement in MGCl would have on the 
metallicity estimates we obtained from scaled-solar models. 
The two a-enhanced BaSTI isochrones c orrespond to total 
metallicities of [M/H] = -2.27 and -1.79 (|Pietrinferni et al.1 
I2OO6I ). and therefore have lower iron abundances than the 
scaled-solar models in Fig. [S] The three a-enhanced Dart- 
mouth models, on the other hand, have the same iron 
abundances as the previous scaled-solar models and there- 
fore higher t otal metallicities. Adopting the relation of 
ISalaris et al] l| 19931 ): 

[M/H] = [Fe/H] + log(0.638/<, + 0.362) , (5) 
where fa = ]^q[«/f<=1 ^j^g a-element content, these 
isochrones correspond to total metallicities of [M/H] ~ —2.2, 
— 1.9 and —1.6, respectively. From Fig. [7] it is clear that 
for both the BaSTI and Dartmouth models it is the most 



metal poor isochrone that provides the best fit to the clus- 
ter RGB, implying that if MGCl has an a-enhanced rather 
than solar- like composition it is still a very metal poor object 
with [M/H] « —2.2 or —2.3. Once again, this is in excellent 
agreement with the metallicity estimate that we obtained 
by aligning Galactic globular cluster fiducials. 

For each set of isochrones we again used the offsets 
required to align the best-fitting model with the CMD to 
determine estimates of the MGCl distance modulus and 
foreground extinction. Our results are presented in Table 
[2I while Fig. |8] shows examples of the aligned best-fit scaled- 
solar isochrones from both group. As previously, random 
uncertainties in the derived distance moduli and colour ex- 
cesses are related to the errors associated with aligning the 
isochrones correctly, and we assume these to be better than 
0.02 mag in E{B - V) and ±0.1 mag in p.. 

The derived values of distance modulus and foreground 
reddening listed in Table [2] show good consistency irrespec- 
tive of the adopted isochrone set or a-element abundance. 
More specifically, the calculated distance moduli span the 
range 23.93 - 24.10 with an rms value of /.i = 24.00 ± 0.06, 
while the E{B — V) measurements span the range 0.15 — 0.19 
with an rms value of 0.17 ± 0.02. These results are entirely 
consistent with those we derived by fitting Galactic globular 
cluster fiducials to our CMDs. 
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Table 2. Results from fitting theoretical isochrones to the CMD of MGCl. 



Isochrone Set 


Filters 


[Fe/H] 


[a/Fc] 


[M/H] 




E{B - V) 


BaSTI 


(9. i) 


-2.27 


0.0 


-2.27 


23.97 


0.18 


BaSTI 




-2.62 


+0.4 


-2.27 


23.93 


0.19 


BaSTI 


(a, r) 


-2.27 


0.0 


-2.27 


24.10 


0.15 


BaSTI 


(9, r) 


-2.62 


+0.4 


-2.27 


24.03 


0.17 


Dartmouth 




-2.35 


0.0 


-2.35 


23.95 


0.18 


Dartmouth 


(9. i) 


-2.50 


+0.4 


-2.21 


24.01 


0.18 


Dartmouth 


(a, r) 


-2.35 


0.0 


-2.35 


24.00 


0.16 


Dartmouth 


{9. r) 


-2.50 


+0.4 


-2.21 


24.08 


0.15 



3.4 MGCl distance and foreground extinction 

Overall, our best estimate for the distance modulus of 
MGCl is 23.91 ± 0.06 from fitting Galactic globular cluster 
fiducials to our CMD, and 24.00 ± 0.06 from fitting theoret- 
ical isochrones. These correspond to line-of-sight distances 
of 605 kpc and 630 kpc respectively, indicating that MGCl 
lies substantially nearer to us than do the central regions 
of M31, which have a line-of-sight d istance of « 780 kpc 
(/iM3i = 24.47 ± 0.07: iMcConnachie e t al. 2005). This result 
is consistent with the p revious distance estimate for MGCl 
bv lMartin et al.l (|2006l ). 

The distance we have measured for MGCl allows us to 
estimate the three-dimensional (deprojected) distance be- 



tween this cluster and the centre of M31. In doing so, we 
adopt /1M31 = 24.47 ± 0.07 along with a corresponding pro- 
jected distance between the cluster and M31 of Rp — 117 
kpc. The derived deprojected galactocentric radius of MGCl 
is then 205 kpc for a distance modulus of 23.91 and Rgc = 
185 kpc for fi — 24.00. Recalling that the la uncertainties 
on our measured distance moduli are ±0.06 mag, the uncer- 
tainties in these deprojected radii are ±20 kpc. 

Meanwhile, our estimate of the MGCl colour excess is 
E{B-V) = 0.18±0.02 from fitting Galactic globular cluster 
fiducials to our CMD, and E{B - V) ^ 0.17 ± 0.02 from 
fitting theoretical isochrones. These values are a factor of 
two larger than the prediction from the reddening maps of 
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Figure 9. Subaru/Suprime-Cairi (/-baud cut-out of MGCl, spanning 6' X 6', showing the spatially variable nebulosity found in this part 
of the sky. The image has been smoothed slightly with a Gaussian kernel of FWHM ~ 2 pixels (fs; 0.4") and then had an appropriate 
linear scaling applied to show the local excess nebulosity clearly. North is towards the top of the image and east is to the left. 



ISchlegel et al] l|l998h for this region of the sky: E{B — V) = 
0.086. Because this discrepancy exists for all of the different 
Galactic globular cluster fiducials we considered and for all 
of the different isochrone sets we used, and because it occurs 
consistently on both the (g, r) and {g, i) CMDs, it is unlikely 
to be a result of our data reduction or analysis. It also cannot 
be due to the omission of colour terms in determining our 
photometric calibrations, as we discuss below. 

A possible explanation for our large inferred E{B — V) 
value may be seen from examination of wide-field imaging 
of the region surrounding the cluster. Fig. [5] shows a 6' x 6' 
g-band image centered on the cluster, excised from Sub- 
aru Suprime-Cam observations targeting the dwarf galaxy 
And XIII (Martin et al. 2009, in prep.) which lies projected 
nearby to MGCl. The image shows nebulosity varying across 
the region, implying that there is variable extinction present 
on scales much sma ller than the 6' x 6' resolution of the 
ISchlegel" et al.l (|1998D maps. In particular, inspection of Fig. 
[9] reveals that the excess nebulosity happens to be large on 
the cluster, consistent with our derived colour excess. We 
note that a similar pattern is also visible on the MegaCam 
survey imaging of this region from which MGCl was discov- 
ered, and also, to some extent, on our CMOS imaging. This 
consistency indicates that the nebulosity visible in Fig. |9] is 
a real feature of this part of the sky, and not simply due to a 
fiat-fielding error or to scattered light on the Subaru image. 

It is worthwhile considering the effects that the omis- 
sion of colour terms from our photometric calibration (see 
Section 12. 2p might have on our derived distance modulus 
and foreground extinction. These quantities were measured 
by registering fiducial sequences to the cluster HB luminos- 



ity and RGB colour. On the {g,i) CMD, the level region 
of the HB occurs at {g — i) ~ 0.5. Eq. [3] then suggests 
that leaving out colour terms in our photometric calibra- 
tion results in an i-band HB magnitude that is too faint 
by ~ 0.01 mag. In other words, incorporating this correc- 
tion would make the distance modulus shorter by ~ 0.01 
mag. Similarly, the colour of the RGB at the HB level is 
{g — i) ~ 1.0. Again using Eq. [31 at this colour the effect of 
omitting colour terms from our photometric calibration is 
negligible, A{g-i) ^ 0.00. For the RGB above the HB level, 
the effect of omitting colour terms is to produce an RGB 
that is too blue by up to ~ 0.01 mag. Thus, including the 
colour terms would only act to increase our derived colour 
excess. Following the same process for the {g,r) CMD, we 
obtain very similar results. Overall, it is then clear that (i) 
including colour terms in our photometric calibration would 
not have resulted in our deriving a larger distance modu- 
lus or smaller foreground extinction for MGCl; and (ii) the 
expected corrections are, in any case, negligible in size. 



4 CLUSTER STRUCTURE AND LUMINOSITY 

We determined the radial density profile of MGCl in each of 
our three photometric passbands via a combination of sur- 
face photometry in its unresolved central regions, and star 
counts at larger radii. We obtained surface photometry by 
using the phot task in iraf to measure the fiux through 
circular apertures of increasing size, and then converting 
these into measurements of the surface brightness in con- 
centric circular annuli. An estimate of the sky background 
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was determined from several blank regions near the cluster, 
and subtracted from each annulus. We desired to extend 
the integrated-light photometry as far from the cluster cen- 
tre as possible to maximize overlap with star-count profiles 
(see below); in practice we were limited to a maximum ra- 
dius Tp = 35", beyond which the flux contribution from ob- 
jects not belonging to the cluster became non-negligible (as 
determined from the CMD). Even within this radius there 
were several significant background galaxies evident on the 
images (see Fig. [l]); these were masked prior to our mea- 
surements. The surface photometry was transformed onto 
the SDSS photometric system via Eq. [5] 

We employed star counts to extend the radial profile 
beyond 35". To minimize contamination from non-cluster 
members we only counted stars falling in several empirically- 
determined regions on both the [g, i) and {g, r) CMDs. On 
each CMD we used stellar detections within 40" of the centre 
of MGCl to define three boxes preferentially inhabited by 
cluster members - one region enclosing the HB, and one 
each for the lower and upper RGB. Because we suspected 
variable foreground extinction on scales of a few arcminutes 
or less (see Section [3]4]), we ensured these CMD regions were 
sufficiently broad to accommodate changes in E{B—V) of up 
to ±0.1 mag away from the cluster centre. Example selection 
boxes for the {g, i) CMD can be seen in Fig. 1101 

Our star-count annuli extended from the innermost re- 
gions of MGCl out to the edges of our stacked images at 
Vp — 150". To ensure consistency with the inner surface 
photometry, in each individual annulus we converted the 
summed fiux of the counted stars into a surface brightness. 
Because this process did not include unresolved diffuse clus- 
ter light, there was a clear offset between our star-count 
profile and that from the inner surface photometry. To join 
the two profiles we considered an overlap region, extending 
inwards from 35" to a radius where we judged that incom- 
pleteness due to crowding was having a significant effect on 
the star-count profile. This was clearly visible as a point 
near rp « 13" where the shapes of the two profiles began 
to diverge - that from the star counts flattening and falling 
below that from the surface photometry. 

Although our CMD selection boxes greatly reduced con- 
tributions to the star-count profile from non-cluster objects, 
there was still undoubtably sufficient contamination that 
it was necessary for us to correct for this. We originally 
planned to use the outermost regions of our CMOS imaging 
to define the background level; however, to our surprise the 
star-count profile was still clearly falling steeply as a func- 
tion of radius to the very edges of the field - indicative of a 
signiflcant contribution from cluster members even at large 
radii. Evidence of this may be seen in Fig. 1111 the left panel 
of which shows a CMD for all stellar objects at radii beyond 
40". Objects falling in the selection boxes on both the ((/,i) 
and (g,r) CMDs are marked with large solid circles. Al- 
though there is certainly contamination in the (g, i) CMD 
selection boxes from non-cluster members, there is also a 
signiflcant number of stars falling very close to where the 
main cluster sequences are expected. In particular a clear 
overdensity of stars is visible at the blue end of the HB, as 
well as around the join between the lower and upper RGB. 
Many of these objects fall at radii even beyond 80". A plot of 
their spatial positions (right panel of Fig. Illjl reveals them 
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Figure 10. All detected stars within 40" of the centre of MGCl 
plotted on a (5,1) CMD. These sequences define the three se- 
lection boxes (marked) that we used to identify probable cluster 
members at larger radii. 

to be fairly evenly distributed about the cluster ~ no striking 
azimuthal asymmetry is evident. 

To obtain an estimate of the background level and trace 
the outermost regions of the cluster profile we utilised two 
supplementary wide-field data sets - photometry from the 
CFHT/Meg aCam survey imaging in which MGC l was dis- 
covered fsee lMartin et al.ll20od : llbata et al.ll2007l '). and pho- 
tometry from the Subaru Suprime-Cam observations noted 
in Section [33] (see Martin et al. 2009, in prep.). The CFHT 
observations are in the MegaCam g and i passband^ and 
are roughly 1.5 — 2 mag shallower than our GMOS obser- 
vations. The Subaru observations are intermediate in both 
depth and image quality between our GMOS observations 
and the CFHT imaging. Both the Subaru/Suprime-Cam and 
CFHT/MegaCam photometric data sets were obtained as 
end-products from ve rsions of the CASU reduction pipeline 
(I Irwin fc Lewisll200ll ). Both cover regions of sky to radii at 
least ~ 2000" from the centre of MGCl. 

To utilize these data sets we followed a similar proce- 
dure to that for constructing our GMOS star-count profile. 
We first defined appropriate selection boxes on the respec- 
tive (3, i) CMDs by using stars within 40" of the cluster 
centre, although we were careful to (i) exclude the most 
crowded inner regions in both sets where the photometry 
is degraded, and (ii) not include objects too close to the 
faint limit in the Subaru photometry, which is spatially vari- 
able over the full field of view. We then used these selection 
boxes to perform star counts in concentric circular annuli 
extending to 400" from the cluster centre. We determined 
the background level by considering a wide annulus exte- 
rior to this, spanning rp = 400 — 2000". We joined the two 
profiles to our composite GMOS profile by considering the 

^ Which are similar, but not identical, to the SDSS g and i filters. 
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Figure 11. All stellar objects in our GMOS field outwith 40" from the centre of MGCl. Left panel: Plotted on a (g, i) CMD along with 
the selection boxes from Fig. 1101 Right panel: Spatial positions on the GMOS field of view. In each panel, stars which fell within selection 
boxes on both the (g, i) and (g, r) CMDs are marked with solid black circles; all other stellar objects are marked with open circles. In 
the left panel there are clear overdensities of stars visible at the blue end of the HB and halfway up the RGB, indicating the presence of 
cluster members at large radii. The right panel shows that these members extend to the edge of the field of view (marked with dotted 
lines) and appear to be fairly evenly distributed about the cluster centre (marked with a cross). 



overlap in the radial range 20 — 150". Note that the Mega- 
Cam photometry is not strictly on the same system as the 
GMOS and Subaru photometry; however assuming there is 
no strong colour gradient in the outer regions of the MGCl 
profile, the difference is, to a good approximation, a small 
constant offset in magnitudfl. 

The final p-band background-subtracted radial bright- 
ness profile for MGCl may be seen in Fig. 1121 The back- 
ground levels measured from the MegaCam and Subaru pho- 
tometry agree to within 0.1 mag and are marked by the hor- 
izontal dotted lines. We determined uncertainties in the cen- 
tral surface photometry by splitting a given annulus into up 
to 8 segments and considering the variance in flux between 
them; the resulting error-bars are smaller than the sizes of 
the plotted points. Uncertainties in the star count profiles 
were calculated using a combination of Poisson statistics 
and the uncertainty in the background subtraction. To avoid 
congestion in the Figure we have marked error-bars on the 
Subaru annuli only, however these are indicative of the typ- 
ical uncertainty at a given radius. The vertical dashed line 
marks the FWHM of the GMOS g-band PSF and is thus 
the approximate radius interior to which atmospheric see- 
ing strongly affects the cluster profile. 

The GMOS measurements alone (solid black circles 
and magenta triangles) show that the cluster can be reli- 
ably traced to a radius rp = 150" ~ 450 pc (assuming 
pi = 23.95), at which point the surface brightness is com- 
parable to the background contribution. Points exterior to 
this are strongly sensitive to the adopted background level, 
but, taken at face value, the profile apparently extends to 



\http://www2. cade- ccda.hia-iha.nrc-cnrc.gc.ca/megapipe/ docs/ filters. html\ 



~ 250 - 300" « 750 - 900 pc without showing a clear tidal 
cut-off. 

We attempted to fit, by minimization, two sets of 
commonl y-used mode ls to the MGCl profile. The first were 
empirical iKiiiS l| 19621 ) models of the form: 

2 

1 1 ' 



E(r, 



So 



(6) 



where rt is the tidal cut-off radius, rc is the core radius, 
and So is a scaling parameter. As long as one is not dealing 
with a low-concentration cluster (ie, provided rt ^ rc), Eo is 
approximately equal to the central surface brightness, and 
Tc is roughly the radius where the surface brightness has 
fallen to half its central value. 

We also fit models of the form described by 
lElson. Fall fc FreemanI l|l987l . hereafter EFF): 



E(.p)=Eo(l + ^ 



(7) 



which are similar to the lKin j (|l962l ) models in their central 
regions but have no outer tidal cut-off, instead approaching 
a power-law fall-off with exponent —7 at large radius. Here, 
the scale radius a is related to the King core radius rc by: 



aV2T_l (8) 

provided, once again, that rt ^ r^- 

The results of our model fitting are displayed in Fig. 
1121 The left panel shows that while a King profile fits the 
inner regions of MGCl well, in the outer parts these models 
do not provide a good description of the observed shape of 
the profile, even if a wide variety of tidal radii are invoked 
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Figure 12. Radial background-subtracted SDSS g-band surface brightness profile for MGCl. Tiie solid black circles indicate integrated- 
light photometry in the unresolved central regions from our GMOS imaging, while the solid magenta triangles denote star-count mea- 
surements at larger radii from the same observations. At larger radii still, star counts from our CFHT/MegaCam survey imaging are 
indicated by red stars, while star counts from Subaru/Suprime-Cam are marked by open blue circles. Uncertainties are plotted on the 
Subaru points; at given radius these are typical for all the star-count measurements. Background levels derived from the MegaCam and 
Suprime-Cam observations are marked with horizontal dotted lines; these agree to better than 0.1 mag. The vertical dashed line ind icates 
the 0. 5" FWHM of the g-band PSF; at smaller radii atmospheric seeing strongly affects the profile. The left panel shows three iKind 
lll962l) models with = 0.65" and n = 35", 200" and 1000". The right panel shows our best-fitting EFF model, which has a = 0.8" 
and power-law fall-off 7 = —2.5. There is a break to a steeper power-law fall-off with 7 = —3.5 (long-dashed line) around rp ~ 35". 



(the three examples plotted have rc = 0.65" and Vt — 35", 
200" and 1000"). The EFF family of models does somewhat 
better, although a single set of parameters cannot describe 
the full radial extent of the observations. The cluster appar- 
ently closely follows a power-law fall-off in its outer regions, 
but there is a break around rp ~ 35" from a profile with 
exponent 7 « —2.5 to a steeper fall-off with 7 « —3.5. It is 
notable that the break occurs roughly where our composite 
profile transitions from surface photometry to star counts, 
and a stellar luminosity function that steepens with radius 
could thus potentially be responsible. However, the fact that 
the break can apparently be seen when considering only the 
GMOS star counts alone (magenta triangles) argues against 
this interpretation. 

Our best-fit models allow us to precisely estimate the 
cluster core radius: the King profiles have rc — 0.65" « 2.0 
pc, while the EFF model has a — 0.80" which corresponds 
to Tc = 0.69" « 2.1 pc (assuming 7 = —2.5, which is 
most appropriate for the inner part of the observed profile). 
Both sets of models converge to a central surface bright- 
ness Sg(0) — 17.85 mag/arcsec^. Our measured values for 
rc are close to the FWHM of the PSF for the g-band ob- 
servations, and could therefore conceivably be smaller than 
we have measured here. Such a scenario would also result 
in a larger central surface brightness than we have deter- 
mined. Higher resolution observations of MGCl will be re- 
quired to test this possibility. For now we simply note that 
many globular clusters with smaller rc and higher central 



surface brightness than we have measured for MGCl are 
known both in the inner and outer regions of M3 1 (e.g., 
iBarmbvT Holland fc Huchrall200 j : iTanvir et al.ll2009l ). 

Our observed profile also allows an estimation of the 
cluster luminosity. Integrating the measured surface bright- 
ness out to the edge of the GMOS field at rp — 150", 
and subtracting the Subaru/MegaCam background level re- 
sults in a total apparent luminosity rrig = 15.743- Assuming 
fi = 23.95±0.06 and E{B-V) = 0.18±0.02 leads to an abso- 
lute g-band magnitude Mg — — 8.9±0.2. We can convert this 
to a standard V-band magnitude by adopting Mv,q ~ +4.82 
and Mg,Q = -1-5.12 (see the SDSS DR7 website as spec- 
ified in Section [22J. This results in Mv = -9.2 ± 0.2, 
rendering MGCl a rather luminous globular cluster with 
itot ~ 4.05 X 10^1/0. Our result is significantly mor e lumi- 
nous than the value derived by iMartin et al.l (|2006l ) : how- 
ever the difference can be entirely attributed to our larger 
assumed foreground extinction, as well as the fact that we 
have traced the cluster profile to a much larger radius. 

Repeating this exercise for the r- and i-bands leads to 
apparent integrated colours of (g — r) = 0.64 and {g — i) = 
1.03. Correcting for foreground extinction and transform- 
ing to standard magnitudes using the relations on the SDSS 
DR7 webpage results in (1/ - J)o = 0.89 for MGCl, which is 



^ It is unnecessary to integrate to larger radii, where the 
background-subtracted profile is quite uncertain, as this adds at 
most a few hundredths of a magnitude to the total luminosity. 
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completely consistent with the integrated colours of other re - 
mote globular clusters in both M31 (e g- jHuxor et al.ll2009l ) 
and the Milky Way fe.g.. lHarrij |l996'). 

Finally, returning to our g-hand profile it is straight- 
forward to determine the half-light radius rh = 2.5" ~ 7.5 
pc . This is larger than the value of 2.3 ± 0.2 pc measured 
bv lMartin et al.l (|2006l ): however much of this discrepancy is 
due to our larger m easured integr ated luminosity - adopting 
that determined bv lMartin et aLi (2006) yields an rh smaller 
by at least a factor of two. A summary of our structural and 
luminosity measurements for MGCl is provided in Table [S] 



5 DISCUSSION 

The measurements we present above suggest that MGCl lies 
at a distance of 200 ± 20 kpc from the centre of M31. This 
renders it by far the most isolated globular cluster known 
in the Local Group - it resides at a much greater galac- 
tocentric radius than do the most remote globular clusters 
in the Milky Way (Palomar 4 at R^c ~ 110 kpc and AM 
1 at Rgc ~ 120 kpc) and the next most distant globular 
clusters known in M31 (e.g., GC10/H 27 at Rp ~ 100 kpc, 
iMackev et"aIll2007l : lHuxor et al.ll2008l l. Indeed, our derived 
galactocentric radius for MGCl is likely to represent a sig- 
nificant fraction of the virial radius of M31. Even so, the 
cluster's radial velocity of ~ —50 km s~^ relative to M31 
is consistent with that expected for a halo member at its 
observed galactocentric radius, and lies wi thin the proba- 
ble M 31 escape velocity (see e.g.. Fig. 4 in lChapman et ahl 
I2OO7I ). There is no implication, therefore, that MGCl is any- 
thing so exotic as an unbound or intergalactic globular clus- 
ter. Instead, our observations suggest that the globular clus- 
ter populations of large spiral galaxies may be considerably 
more spatially extended than has previously been appreci- 
ated. 

In this regard it is interesting that although distinct 
substructures have been identified in the M31 halo to 
projected radii beyond 100 kpc (e.g., Ilbata et"aD l2007l : 
iMcConnachie et al.l |2009| ). MGCl does not lie near any 
known tidal stream or halo overdensity. Rather, it falls 
within an area that is apparently com pletely free of such 
structures - indeed it is this region that Ilbata et all (|2007l ) 
use to define their "clean" M31 outer-halo stellar selection. 
It is also worthwhile noting that while MGCl lies within a 
~ 2° X 2° area that contains three satellite dwarf galaxies of 
M31 (And XI, XII and XIII), there is no evidence that any 
of these four objects is associated with any other. On the 
contrary, the three dwarfs appar ently lie much close r to the 
M31 distance than does MGCl (jMartin et al.ll2006l . 2009 in 
prep.), w hile all four exhibit qu ite disparate radial velocities 
(see e.g.. IChapman et al.ll2007l '). 

Our measurements have also shown that MGCl is a 
very luminous globular cluster. Indeed, with My ~ —9.2 
it would fit comfortably within the top ten most luminous 
globular clusters in the Milky Way. The fact that it also lies 
at a very large galacto centric radius reinforces the idea, first 
put forward by Mackev et al.l (|2007t ) (see also iGalleti et al.l 
l2007l : lHuxor et al.ll2009l ). that M31 possesses a population of 
very remote, compact, luminous globular clusters not seen 
in the Milky Way. A total of seven globular clusters are 
known in the very outer reaches of the Galactic halo, beyond 



Table 3. Structural and luminosity measurements for MGCl. 



Parameter 


Value" 


Sg(0), Sr(0), S,(0) 


17.85™, 17.17™, 16.83' 


Sv,o(0) 


16.95™ 


Tc (King) 


0.65" (f5! 2.0 pc) 


n (King) 


Indeterminate 


rc (EPF) 


0.69" (^ 2.1 pc) 


7 (EPF) 


-2.5 ^ -3.5 


rh 


2.5" (» 7.5 pc) 


rrig 


15.74™ 


Mg 


-8.9™ 


Mv 


-9.2™ 


Ltot 


4.05 X W^Lq 


(g - r), (g - i) 


0.64™, 1.03™ 


(V - I)o 


0.89™ 



" A distance modulus of fi = 24.0 and foreground extinction 
E(B — V) = 0.17 have been used where appropriate. 

Rgc = 40 kpc. Of these, six are subluminous (with —6.0 < 
Mv < -4.7) and rather diffuse (11 < r^, < 25 pc). The sev- 
enth cluster is NGC 2419, which is one of the most luminous 
in the Galactic system and which has a number of unusual 
properties that has led to suggest ions that it may be the re- 
mains of an accreted galaxy (e.g.. Ivan den Bergh fc MackevI 
2004 ). In contrast, a search of the RBC reveals 34 confirmed 
clusters in M31 with Tip ^ 40 x (f ) = 31.4 kpc0, of which 25 
are compact and 9 are of the extended type first reported by 
iHuxor et al.l (|2005l ). Assuming /im3i ~ 24.47 and a typical 
E{B — V) ~ 0.1, 12 of the 25 remote compact clusters (in- 
cluding MGCl) have luminosities brighter than My ~ —8.0, 
while another 5 have luminosities lying between this value 
and the peak of the GCLF at My ^ -7.5. 

Given that all but three of the 34 remote M31 c luster s 
come from the survey area described in iHuxor et al.l (|2008h . 
covering one quadrant of the M31 halo, this clear disparity 
between the remote M31 and Milky Way globular cluster 
systems cannot be due simply to the global difference in 
population size (which would account for, at most, a factor 
~ 3). The origin of the differences between the Galactic and 
M31 outer halo globular clusters is not known. We speculate 
that it may be related to the apparently more vigorous ac- 
cret ion and merg er history that M31 has experienced (e.g., 
Iba ta et al.ll2007l ): however far more detailed study of re- 
mote M31 members will be required to shed light on this 
possibility. 

MGCl offers a useful probe of the outermost regions 
of the M31 halo, wher e only a very few stars have previ- 
ously been m easured jKalirai et al. I [2006; Chap man et al.l 
[iood . I2OO8I : iKoch et al.l |2008[) . The metal abundance we 
derive for MGCl is low, [Fe/H] « -2.3. This is in con- 
trast to results from rec ent high-resolution spectroscopy by 
lAIves-Brito et al.l (j^OOg") who found [Fe/H] = -1.37 ± 0.15. 
The origin of this discrepancy is unclear; however we note 
that it also exists for the other two clusters in their sam- 
ple: GC5 and GCIO, which are the next two most re- 



^ Here the factor ^ scales the three-dimensional radius Rgc = 40 
kpc to an average projected radius. 
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lAlves-Brito et all (|2009l ) derive [Fe/H] = -1.33 ± 0.12 and 
[Fe/H] = —1.73 ±0.20, respectively, for these objects, while 
high-quality CMDs s uRRest [Fe/H] = - 1.84 ± 0.15 and 
[Fe/H] = -2.14±0.15 jMackev et al.ll2007l ). Until these dis- 
agreements are resolved, it is too early to conclude that the 
mean globular cluster metallicity at large radii in M31 flat- 
tens off at [Fe/H] ~ —1.6 as suggested bv I Alves-Brito et al.l 
(|200&y Taking the CMD results at face value would instead 
indicate that [Fe/H] < —2.0 for clusters beyond i?p ~ 80 
kpc. 

The situation for M31 halo stars at radii beyond 
~ 70 k pc is, at present , conf using (see e.g., the discus- 



sion in iRichardson et al 



20091). Spectro s copic studi es by 



iKalirai et ahl (|2006l ) and IChapman erall l|2006l . |2008| ) sug- 
gest ( [Fe/H] ) ~ -1.3 in the range 60 < i?p < 160 kpc, not 
inconsistent with the peak of the (photometric) metallicity 
distribution function (MD F) for outer-halo n on-stream stel- 
lar populations derived bv llbata et al] l|2007l ). On the other 
hand, after reanalysing of much of t he sp ectroscopic data 
from the above studies, iKoch et all (12008') found that the 
mean metal abundance of M31 halo stars may fall below —2 
beyond -Rp ~ 80 kpc, and rea ch as low as ([Fe/H] ) ~ —2.5 
at ilp - 160 kpc. The MDF of llbata et all (|2007l ') appears to 
exhibit a metal-poor tail and our results for MGCl suggest 
that there must be at least some very metal-poor compo- 
nent to the extremely remote outer halo of M31. However 
on the whole it seems that the most distant globular clusters 
in M31 are more metal-poor than the bulk of the field halo 
stars at comparable radii. 

Another remarkable feature of MGCl is its extremely 
extended outer structure. Cluster members are clearly vis- 
ible at large radii, and these are apparently fairly homo- 
geneously distributed azimuthally. We have traced the ra- 
dial profile of MGCl to at least ~ 450 pc and possibly 
as far as ^ 900 pc, rendering it by far the most extended 
globular cluster hitherto studied. The most extended clus- 
ter in the Milky Way is NGC 2419, with n « 200 pc; the 
vast majority of the G alactic population have rt < 100 pc 
fe.g.. lMackev &: van den Bcr gh 2005). The unique structure 
of MGCl likely reflects its unmatched remote location, away 
from any strong tidal forces which might act to impose a 
tight limit on its radial extent. The expected tidal radius of 
a globular cluster of mass M^i orbiting a galaxy of mass Mg 
at a galactocentric radius i?gc is given by: 



V 3Mg 



(9) 



Adopting our measured parameters for MGCl (-Rgc = 200 
kpc and Ltot = 4.05 x IQ^Lq) along with a typical V-band 
mass-to- light ratio for a globular cluster M/L « 2 and the 
meiss of M31 from |e vans fc Wilkinsoiil (|2000l ). Mg ^ 1.2 x 
IO^^Mq, we find the expected tidal radius for MGCl is rt ~ 
1200 pc. It is therefore not surprising we have not identified 
a tidal cut-off in its radial profile. 

Note that we cannot be certain the outermost cluster 
members are bound to the cluster, since any tidal distortions 
would not necessarily lie in the plane of the sky. However, 
the lack of any clear spatially asymmetric distribution of 
these outer members together with the calculation above 



showing a very large expected rt suggests more of a halo- 
type structure than strong tidal tails. 

Assuming a bound halo of stars around the cluster, the 
observed power-law fall-off of the profile is not unexpected. 
iBaumgardt. Hut fc Heggi^ (|2002l ') used A^'-body simulations 
to study the long-term evolution of isolated star clusters. 
Their models show that an isolated cluster tends to build 
up a surrounding halo of stars that have been scattered onto 
radial orbits by two- or three-body encounters in the inner 
regions of the cluster. The resulting density profile follows 
a power-law fall-off as p ~ r~^ '^ at large radii. This is very 
similar to the results fr om, for example, the Monte Carlo 
models of ISpitzeil (|l987t ) which predict p ~ r~^'^ . In projec- 
tion we thus expect a power-law fall-off with an exponent 
somewhere between —2.5 < 7 < —2.3, which is entirely con- 
sistent with what we observe for the inner profile of MGCl. 

In both sets of models the core-halo structure takes 
many median relaxation times to establish over radial ranges 
of up to ~ lOOrh ~ 750 pc for MGCl. It is therefore per- 
haps not surprising that we observe a break to a steeper 
fall-off in the outer regions of MGCl. This may reflect the 
ongoing initial build-up of the cluster's halo, or possibly the 
re-establishment of this halo after the last pericentre pas- 
sage of MGCl which could have brought it much closer to 
M31 than it is at present. Whatever the case, given that 
the median relaxation time of MGCl is likely to be several 
Gyr (s ee, for examp le, the values listed for Galactic globu- 
lars bv lHarrisllTogg ). the existence of a core-halo structure 
at all suggests that MGCl has been evolving in isolation for 
a very considerable period of time. 



6 CONCLUSIONS 

In summary, we have presented resolved photometric mea- 
surements of MGCl from deep Gemini/GMOS imaging of 
this remote M31 cluster. Our resulting CMD displays fea- 
tures consistent with those of an ancient stellar popula- 
tion, including a long horizontal branch populated across 
the instability strip. By fitting isochrones calculated by sev- 
eral different groups, as well as fiducial sequences observed 
for Galactic globular clusters, we determine MGCl to be a 
metal-poor object with [M/H] ~ —2.3. This is in disagree- 
ment wit h a recent high re s olutio n spectroscopic measure- 
ment by lAlves-Brito et ahl (|2009l ) who find a significantly 
higher metal abundance; the origin of this discrepancy is 
not yet clear but it also exists for several other remote M31 
globular clusters. 

The best-fitting stellar evolution models and globu- 
lar cluster fiducials suggest that MGCl lies significantly 
closer to us than M31. Our preferred distance modulus is 
= 23.95 ±0.06 (corresponding to a line-of-sight distance of 
~ 620 kpc), compared with the generally accepted value for 
M31 of PM31 = 24.47 780 kpc). Combined with the large 
projected separation between MGCl and M31 {Rp = 117 
kpc), our distance measurement implies that MGCl lies in 
the extreme outer regions of the M31 halo with a galacto- 
centric radius 7?gc = 200 ± 20 kpc. This renders it the most 
isolated known globular cluster in the Local Group by some 
considerable m argin. Even so, the radial velocity measured 
for this object (|Galleti et all l2007l : lAlves-Brito et al.ll2009l ) 
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is within the M31 escape velocity and consistent with that 
expected for an M31 halo member at the radius of MGCl. 

By measuring the radial brightness profile of MGCl we 
have established that this cluster is centrally concentrated 
with a core radius rc ~ 2 pc. It is also rather luminous with 
M y = —9.2. Th i s rein forces the notion, first put forward 
bv lMackev et al.l (|2007D . that M31 possesses a population of 
compact, luminous globular clusters in its outer halo that 
does not have a counterpart in the Milky Way. 

Perhaps unsurprisingly given its extremely remote lo- 
cation, MGCl does not show any evidence for a tidal limit 
in its outer regions, and King models do not provide a good 
description of the shape of the brightness profile at large 
radii. Instead we observe a power-law fall-off with exponent 
7 — —2.5, breaking to 7 = —3.5 in the outermost parts. 
This core-halo structure, which we see extending to radii of 
at least 450 pc and possibly as far as ~ 900 pc, is broadly in 
line with expectations derived from numerical modelling of 
isolated globular clusters. In this regard MGCl is unique in 
our experience - the most radially extended globular cluster 
in the Milky Way is NGC 2419 which shows a tidal cut- 
off at rt ~ 200 pc. Further detailed observational study of 
MGCl may therefore be extremely useful for advancing our 
understanding of globular cluster evolution. 
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